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Abstract 

In-medium ciiiral perturbation theory is used to calculate the density dependence of 
the quark condensate {qq)- The corrections beyond the linear density approximation 
' are obtained by differentiating the interaction contributions to the energy per particle 

■ of isospin-symmetric nuclear matter with respect to the pion mass. Our calculation treats 
systematically the effects from one-pion exchange (with m,r-dependent vertex corrections), 
iterated Ivr-exchange, and irreducible 27r-exchange including intermediate A(1232)-isobar 
excitations, with Pauli-blocking corrections up to three-loop order. We find a strong and 

. non-linear dependence of the "dropping" in-medium condensate on the actual value of the 

^— ^ I pion (or light quark) mass. In the chiral limit, = 0, chiral restoration appears to be 

■ reached already at about 1.5 times normal nuclear matter density. By contrast, for the 
physical pion mass, = 135 MeV, the in-medium condensate stabilizes at about 60% of 
its vacuum value above that same density. Effects from 27r-exchange with virtual A(1232)- 
isobar excitations turn out to be crucial in generating such pronounced deviations from 

Q . the linear density approximation above po- The hindered tendency towards chiral symme- 

5 I try restoration provides a justification for using pions and nucleons as effective low-energy 

degrees of freedom at least up to twice nuclear matter density. 



> ; PACS: 12.38.Bx, 21.65.+f 

^ ■ Keywords: In-medium chiral condensate, long-range correlations in nuclear matter from one- 
and two-pion exchange. 
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1 Introduction and framework 



The quark condensate |(0|gg|0)| is an order parameter of spontaneous chiral symmetry breaking 
^ I in QCD. With increasing temperature the quark condensate decreases (or "melts"). For low 
?H I temperatures this effect can be systematically calculated in chiral perturbation theory. At three- 
loop order [1] the estimate ~ 190 MeV for the critical temperature, where chiral symmetry will 
eventually be restored, has been found. This extrapolated value of Tc is remarkably consistent 
with Tc = (192 ± 8) MeV [2] obtained in numerical simulations of full QCD on the lattice, 
although this result is still under debate. It has subsequently been criticized [3] with respect to 
the reliability of the continuum extrapolation performed in ref. [2]. In fact the QCD chiral "phase 
transition" is merely a non-singular cross-over, with the transition temperature (deduced from 
the peak of the chiral susceptibility) lying in the broad range Tc = (160 ±29) MeV according to 
ref. [3]. 

The chiral condensate \{qq)\ drops also with increasing baryon density. Presently, it is 
not feasible to study this phenomenon rigorously in lattice simulations of QCD due to the 
problems arising from the complex-valued Euclidian Fermion determinant at non-zero baryon 
chemical potential. As an alternative, the density dependence of the quark condensate {qq){p) 
can be extracted by exploiting the Feynman-Hellmann theorem applied to the chiral symmetry 
breaking quark mass term rriq qq in the QCD Hamiltonian. The leading linear term in the 
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nucleoli density p is then readily derived by differentiating the energy density of a nucleonic 
Fermi gas, pM^ + 0{p^^^), with respect to the light quark mass rriq. This introduces the 
nucleon sigma-term gn — {N\mqqq\N) = niqdMM /druq = (45 ± 8) MeV [4] as the driving 
term for the density evolution of the chiral condensate. Following this simple linear density 
approximation one would naively estimate that chiral symmetry gets restored at (2.5 — 3)po) 
with po = 0.16 fm~^ the nuclear matter saturation density. 

Corrections beyond the linear density approximation arise from the nucleon-nucleon interac- 
tions which transform the nucleonic Fermi gas into a nuclear Fermi liquid. These corrections have 
been studied in one-boson exchange models of the NN-interaction combined with the relativistic 
Dirac-Brucckncr approach to nuclear matter [5, 6]. Knowledge of the quark mass derivatives of 
the various meson masses and coupling constants is required in order to quantify the interaction 
effects on the in-medium condensate. For the vector and scalar bosons it has been assumed that 
their sigma-terms scale linearly with that of the nucleon, i.e. niq dmv,s / drUg — (TNmv^s/^N- 
In the further development it has been demonstrated in ref.[7] that the higher order corrections 
(in density) depend sensitively on the interpretation of the isoscalar scalar "a" -boson (which is 
responsible for the central NN-attraction in one-boson exchange models) and its substructure. 
Within modest variations of an unknown parameter < Cs < 1 both an accelerated and a 
hindered tendency towards chiral restoration are possible. The way out of this dilemma is to 
replace the fictitious "cr" -boson exchange by realistic two-pion exchange processes. 

Because of the Goldstone boson nature of the pion with its characteristic mass relation, 
~ iriq, the pion-exchange dynamics in nuclear matter plays a particularly important role 
for the in-medium quark condensate. A first step in this direction was made in ref. [8] where 
realistic saturation of nuclear matter could be obtained from the iteration of l7r-exchange plus 
a short-range NN-contact interaction to second order. It was found that the pionic interaction 
effects (with well-known quark mass derivative) counteract the reduction of the condensate 
from the leading linear term in density. Actually, if one restricts oneself to the density region 
P < Po = 0.16 fm~'^, then all existing calculations agree that the deviations from the linear 
density approximation are relatively small and practically get masked by the uncertainty of the 
empirical nucleon sigma-term, gn — (45 ± 8) MeV. 

The chiral approach to nuclear matter has been extended and improved in rcfs.[9, 10] by 
systematically including effects from irreducible 27r-exchange together with excitations of vir- 
tual A(1232)-isobars. The physical motivation for such an extension is threefold. First, the 
spin-isospin-3/2 A(1232)-resonance is the most prominent feature of low-energy TrA/'-scattering. 
Secondly, it is well known that the 27r-exchange between nucleons with excitations of virtual A- 
isobars generates the medium- and long-range components of the isoscalar central NN-attraction 
[11], which is simulated by the scalar "o""-boson in the one-boson exchange models. Thirdly, 
the delta-nucleon mass-splitting A = 293 MeV is of the same size as the Fermi momentum 
kfo — 263 MeV ~ 2mT^ at nuclear matter saturation density. Therefore pions and A-isobars 
should both be treated as explicit degrees of freedom in the nuclear many-body problem. It has 
been found in ref. [10] that the inclusion of the chiral ttNA dynamics significantly improves e.g. 
the momentum-dependence of the (real) single-particle potential U {p, kf) and the isospin prop- 
erties (as revealed by the density-dependent asymmetry energy A{kf) and the neutron matter 
equation of state). Moreover, it guarantees spin-stability of nuclear matter [12]. 

The purpose of the present paper is to investigate the in-medium chiral condensate (beyond 
the linear density approximation) in this extended and improved framework for nuclear matter 
where interaction effects are (almost) exclusively given by one- and two-pion exchange according 
to the rules of chiral symmetry. Since the pion mass ttIt^ appears as an explicit variable in our 
calculation we can also study how the "dropping" in-medium condensate evolves from the chiral 
limit, ttItt = 0, to the real world with its fixed amount of explicit chiral symmetry breaking, 
ttItt = 135 MeV. We find that the in-medium condensate behaves very differently in both cases. 



with drastic consequences for nuclear physics in the chiral hmit. 

Our starting point is the Feynman-Hellmann theorem which relates the in-medium quark 
condensate {qq){p) to the quark mass derivative of the energy density of isospin-symmetric spin- 
saturated nuclear matter. Using the Gell-Mann-Oakes-Renner relation m^f^ — —■mq{0\qq\0) one 
finds for the ratio of the in-medium to the vacuum quark condensate: 

mm fiyK\ wmi^ im%>^"^»> 

with the Fermi momentum kf related to the nucleon density p = 2A;j/37r^ in the usual way. The 
term proportional to = {N\mgqq\N) = rUqdMjq /drUq comes from the noninteracting Fermi 
gas including kinetic energy contributions expanded up to order M^^. Interaction contributions 
beyond the linear density approximation are collected in the function: 

defined as the derivative of the interaction energy per particle E{kf) with respect to m^. We 
mention here that denotes the pion decay constant in the chiral limit and stands for the 
leading linear term in the quark mass expansion of the squared pion mass. With this convention 
the Gell-Mann-Oakes-Renner relation, m^/^ = —mq{0\qq\0), becomes exact and (0|ggf|0) is the 
vacuum condensate in the chiral limit. 



2 Pion mass derivative of the interaction energy 

In this section we present analytical expressions for the contributions to the derivative function 
D{kf) as given by various one- and two-pion exchange diagrams. Taking the m^-derivative is 
a straightforward procedure since we can borrow here heavily from the explicit expressions for 
the diagrammatic contributions to E{kf) written down in our previous works [9, 10]. 

2.1 One-pion exchange Fock term 

We are working to three-loop order in the energy density. At that order one encounters also 
pion-loop corrections to the pion-nucleon vertex. As a consequence, the Ivr-exchange Fock term, 
eq.(6) in ref.[9], needs to be multiplied by the following renormalization factor: 

with Qa the nucleon axial vector coupling constant in the chiral limit. The last two terms, 
depending on the dclta-nucleon mass splitting A = 293 MeV, arise from pion-loop diagrams 
involving the A(1232)-isobar. The low-energy constant 7(A) takes care of the (empirical) 
Goldberger-Treiman discrepancy. It is determined for any choice of the regularization scale 
A by the condition (g'7rAr/-/^Ar)phys — Qa/ f-K, where "phys" denotes physical values. 



Fig. 1: One-pion exchange Fock diagram with pion self energy. Its combinatoric factor is 1/4- 



The Ivr-exchange Fock diagram (with unrenormahzed couphng constant) including the rel- 
ativistic l/M|r-correction leads to the expression: 
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where we have introduced the useful dimensionless variable u = kf/mj^. For reasons of consis- 
tency with the loop expansion the renormalization factor r{mT^) = 1 + 0{m1) must be applied 
only to the static term in eq.(4), and the m^-derivative of r(m7r) must be multiplied only with 
the static Ivr-exchange energy per particle. The necessity for this procedure will become clear 
in subsection 2.6 when discussing the chiral limit —>■ 0. 

The loop diagram with a tadpole on the exchanged pion (see Fig. 1) generates a momentum- 
independent pion selfenergy (i.e. a pion mass shift). Its contribution to the function D{kf) 
reads: 
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with the low-energy constant /3(A) determined by the relation = — 647r^Z3(A) — 2 In(m7r/A) ~ 3 
[13]. Although this contribution is very small it has to be kept for reasons of consistency. 

2.2 Iterated one-pion exchange 

The Hartree diagram from iterated Ivr-exchange with two medium insertions (see Fig. 3 and 
eq.(7) in ref.[9]) leads to the expression: 
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and the corresponding Fock diagram with two medium insertions gives: 
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In our way of organizing the many-body calculation, the Pauli blocking corrections are repre- 
sented by diagrams with three medium insertions. The contribution of the Hartree diagram 



with three medium insertions to the function D{kf) can be reduced to a one-parameter integrah 
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On the other hand, one gets from the Fock diagrams with three medium insertions: 
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where we have introduced the auxihary function: 
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and the abbreviations s — xy ^ \Ju^ — + x^y'^ and t — xz ^ \/v? — + x'^z'^. Note that the 
expressions in cqs.(6-9) carry the large scale enhancement factor M^. It stems from the energy 
denominator of these iterated diagrams which is proportional to the difference of small nucleon 
kinetic energies. 

2.3 Irreducible two-pion exchange 

The irreducible two-pion exchange with only nucleons in the intermediate state leads to the 
following contribution: 
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as obtained by differentiating eq.(14) in ref. [9] with respect to at fixed kf. Wc have ar- 
ranged the M^-terms in the last line such that the low-density expansion of eq.(ll) starts as 
fcj[ln(m7r/A) + 1/2] with no further additive (regularization-scheme dependent) constant to the 
chiral logarithm. Or stated differently, the underlying 27r-exchange interaction at zero momen- 
tum transfer has been restricted to the non-analytical piece proportional to m^ln(m7r/A). The 
dependence of the contribution in eq.(ll) on the regularization scale A will be discussed in 
section 3.2. 



2.4 Two-pion exchange with virtual A-isobar excitation 

We give first the three-body contributions. The Hartree diagram (see Fig. 2 and eq.(5) in ref. [10]) 
leads to the following closed form expression: 
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while the contribution of the Fock diagrams can be represented as a one-parameter integral: 

(13) 

with the two auxiliary functions Gs{x,u) and Gt{x,u) written down explicitly in eqs.(7,8) of 
ref.[10]. As in our previous works we use the value 3/y/2 for the ratio between the nNA- and 
vrA^A^-couphng constants. 

The two-body terms from 27r-exchange with virtual A-excitations fall into two classes: the 
dominant terms scaling reciprocally with the mass splitting A = 293 MeV, and the remaining 
ones with a more complicated A-dependence. The contribution of the dominant two-body terms 
to the function D{kf) can be written again in closed form: 

43 + 161m2 2, 43 2Slu 437w=^l 
5 ln(l-t-M2)H }. (14) 

It has been derived from the following isoscalar central and isovector tensor one-loop NN- 
scattering amplitudes [11]: 
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with q the momentum transfer between the two nucleons. In these expressions we have kept the 
polynomial pieces proportional to odd powers of the pion mass nij^. These non- analytic terms 
in the quark mass rriq are a unique feature of the chiral pion-loop dynamics. For the remaining 
two-body terms we employ the spectral-function representation [10] and differentiate directly 
the imaginary parts of the TrA^A-loop functions with respect to m^. This gives: 
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where the ^f^-terms in the curly bracket stem from box diagrams and the ^f^-independent ones 
from the triangle diagram. The spectral-function representation in eq.(17) involves one subtrac- 
tion of a term linear in density p — 2A;|/37r^. The associated subtraction constant receives also 
contributions from the pion loop diagrams with a non-analytical dependence on the quark mass. 
We reinstore these distinguished pieces from the chiral pion-loop dynamics by the additional 
contribution linear in density: 



which has the property that it vanishes in the chiral limit. 




Fig. 2: Hartree and Fock diagrams related to the chiral inrNN contact vertex proportional to C\ . 
The combinatoric factors of these diagrams are 1/2 and 1, in the order shown. 

2.5 Chiral ttttNN contact vertex proportional to ci 

In chiral perturbation theory the nucleon sigma-term ctat has a specific nonhnear dependence 
on the quark mass niq. The leading linear term comes from the effective Lagrangian C)^^ = 



Ci NNtTx+ = 4cim^iVA^y 1 — 7?^//^ and the nonlinearities arise from pion-loops and higher 
order counterterms. Putting together all (numerically) relevant pieces the ratio ais[/m\ entering 
eq.(l) takes the form: 
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By chiral symmetry, the ci-term in the effective Lagrangian generates also an additional nnNN- 
contact interaction with vertex insertion: —4dcim?^f~'^5°'^. On the one hand it makes up the 
sizeable logarithmic loop correction to (XAr/m^ in eq.(19). On the other hand it gives rise to an 
additional two-pion exchange contribution to the NN-interaction, and moreover it generates a 
long-range three-nucleon force (see diagrams in Fig. 2). The two-body terms in nuclear matter 
lead to the following contribution to the m^-derivative of the energy per particle: 
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which has been derived from the isoscalar central one-loop NN-scattering amplitude: 
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In addition, there are Ci-contributions from the three-body Hartree diagram: 
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and the three-body Fock diagram: 
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with the auxiliary function G(x,m) defined in eq.(lO). 



2.6 Chiral limit 



The m^-derivative of the energy per particle, D{kf), presented in the previous subsections 
takes a particularly simple form in the chiral limit = 0. In that limiting case almost 
all integrals can be solved and the dependence on the Fermi momentum kf becomes simply 
powerlike (with an exponent determined by the mass dimension of the prefactor) . The subscript 
denotes the function D{kf) in the chiral limit 0. The l7r-exchange Fock term (with the 

renormalization factor in the chiral limit, r(0) = 1) gives: 

and the total contribution from iterated l7r-exchange reads: 
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The contribution from irreducible 27r-exchange: 
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has a logarithmic singularity In(m7r/A). It gets exactly canceled by the m^-derivative of the 
renormalization factor r(m^) applied to the static l7r-exchange: 
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This crucial feature instructs us that one must work consistently with the parameters in the 
chiral limit as they are given by the effective chiral Lagrangian. It is mandatory to include 
only those renormalization effects to physical parameters which are actually generated by the 
pion-loops to the order one is working. The constant C represents the additional effect of a NN- 
contact coupling linear in the quark mass m^, and 7 is a left-over from the Goldberger-Treiman 
discrepancy. The two- and three-body terms from 27r-exchange with single A-excitation scaling 
as 1/A read together: 

where the two-body term eq.(14) has contributed 367r^/35 to the numerical factor in brackets. 
For the remaining two-body Fock terms from 27r-exchange with A-excitations the spectral- 
function representation turns into: 



2x + llx^ + 



(1 + ^2) 



2 



-|- (13a; — 2) arctanx 



Ax + Qx-^ , \ ki 

-\ 4 arctana; > $ 



J VxA 



with the auxiliary function: 
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For low densities the contribution in eq.(29) behaves as /cy ln(A;//A). Finally, the total contri- 
bution from the 7r7rA'"A'"-contact vertex proportional to ci reads, in the chiral limit: 

W"-^f?;4). (31) 
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Fig. 3: The nuclear matter saturation curve underlying our calculation of the in-medium quark 
condensate. Apart from the interaction contributions described in refs.[9, 10] and those propor- 
tional to Ci, it includes one single adjusted term linear in density, E{kf)^^^^ = — 7.64GeV~^ A;^. 

3 Results 

This section presents and discusses results for the in-medium chiral condensate as a function of 
baryon density p. Apart from just collecting and computing the series of terms given in section 
2, this includes also a detailed investigation of the pion mass dependence of the quark condensate 
at given density p. It is furthermore necessary to estimate the quark mass dependence of the 
NN-contact term (i.e. the size of the parameter C introduced in eq.(27)) which encodes short- 
distance dynamics not controlled by the underlying chiral effective field theory. At this point we 
can now take very recent computations of the NN-potential from lattice QCD for orientation. 

3.1 Parameter fixing 

First, we have to fix the parameters. The pion decay constant in the chiral limit f^^ is determined 
by the relation: /jr.phys = fni^ + hini-n/'^T^f-nY] = 92.4 MeV. Choosing the central value 1^ = 
4.4 ± 0.2 of ref.[13] one gets /tt = 86.5 MeV. For the nucleon axial vector coupling constant qa 
in the chiral limit we take the value qa — 1.224 as obtained recently via chiral extrapolations 
of lattice data in ref.[14]. A similar analysis [15] gives for the nucleon mass in the chiral limit 
Mat = 882 MeV and for the low-energy constant ci = — 0.93GeV~^ (as central values). For 
the delta-nucleon mass splitting we take the empirical value A = 293 MeV. This is consistent 
to the order in the loop expansion we are working here. The parameter 7 introduced in eq.(3) 

is determined by the relation: {gnN/Mis[)phys = \jT{mT^) gA/ fn- Taking for the left hand side 
13.2/(939 MeV), we deduce 7 = -1.505 at the regularization scale A = = 882 MeV. Or 
stated differently, for our choice of parameters {gA and /tt in the chiral limit) the Goldberger- 
Treiman relation is exact within one percent. We neglect the 1.2%- difference between the 
physical pion mass m7r,phys = nT'^i^ — hinT'Tr/^T^fnT] (with Z3 ~ 3 [13]) and the leading order one. 
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Fig. 4-' The nucleon-nucleon potential in the ^Sq channel from lattice QCD [16] for three different 
pion masses, tUt^ — (380, 529, 731) MeV. 

ttItt, since this difference is mudi smaller than the splitting between charged and neutral pion 
masses. 

With these fixed parameters we obtain (for the physical value of the pion mass m^r — 
135 MeV) the nuclear matter equation of state as shown in Fig. 3. Besides the Itt- and 27r- 
exchange contributions described in refs.[9, 10] and those proportional to Ci, it includes one 
single adjTistcd term linear in density: E(kfY^'^^^ = — 7.64 GeV~^ A;^. We interpret its strength 
parameter to subsume all unresolved short-distance NN-dynamics relevant for nuclear matter 
at low and moderate densities < p < 2.5po = 0.4 fm~^. Its (weak) implicit quark mass 
dependence will be estimated below. The nuclear matter compressibility K — kj^E'^kfo) 
related to the curvature of the saturation curve at its minimum comes out as X = 295 MeV. 
This value lies at the high side of present semi-empirical determinations of K [17]. 

Next, we estimate the parameter C introduced in eq.(27) which represents quark mass 
dependent effects from the short-range NN-interaction. The simple model of a;(782)-meson ex- 
change would give C^'^^ ~ —0.7, a large correction, choosing an a;A^-coupling constant of order 
10 and the constituent quark counting rule dmi^/dmq — 2. However, since the a;-exchange 
model has no clear connection to the short-range NN-dynamics of QCD, wc follow another 
option. Recently, the nucleon-nucleon potential has been studied within lattice QCD [16] us- 
ing the quenched approximation. The short-distance part (r < 0.6 fm) of this potential in 
the ^So channel is shown in Fig. 4 for three different pion masses, = (380, 529, 731) MeV. 
We can identify the volume integrals Iq — 47r jQ°drr'^V{r) over the repulsive cores of these 
potentials with the strength of a contact-coupling in the ^Sq channel. From the three values 
/q = (83.6, 53.3, 21.9)MeVfm"^ we obtain a mean value for the derivative with respect to the 
squared pion mass: dlo/dm^ ~ — 0.17 GeV~^fm^. Via the contribution E^kfY^^^ = lokj/An'^ 
to the energy per particle we estimate the parameter C as C = ATr'^f^dlo/dml ~ —0.05. In 
comparison to the a;-meson exchange this is a rather small number. These considerations raise 
also some doubts concerning the significance of the vector boson exchange phenomenology as a 
valid picture of the short-distance NN-dynamics. Even with a factor 2, to include an equally 
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Fig. 5: The ratio of the in-medium chiral condensate to its vacuum value as a function of the 
nucleon density p for three different values of the pion mass, rrLj^ = (0, 70, 135) MeV. The 
dashed line corresponds to the linear density approximation using the empirical central value 
aN = i5MeV [4]. 

strong contribution from the ^5"! channel [16], the value C ~ —0.1 affects the condensate ratio 
at nuclear matter saturation density po — 0.16 fm~^ (corresponding to kfo — 263 MeV) only at 
the 3 permille level (and 4 times as much at 2po)- 

We can therefore conclude that the short-range NN-dynamics as given by lattice QCD [16] 
has a negligible influence on the in-medium chiral condensate {qq){p)- The deviations from the 
linear density approximation are primarily caused by the long- and intermediate range Itt- and 
27r-exchange dynamics. 

There is some residual dependence on the regularization scale A left over which is not bal- 
anced by the parameters /3(A) and 7(A) (namely from the last term in eq.(ll)).^ Varying 
A between 0.6 GeV and 1.2 GeV changes the condensate ratio at po by 3.5%. Since this is 
much smaller than the effect induced by the uncertainty of the empirical nucleon sigma-term 
aN = (45 ±8) MeV [4] we stay with the "natural" choice of A = = 882 MeV. When inserting 
into eq.(19) it reproduces correctly aj^ — 44.3 MeV for — 135 MeV. 

3.2 In- medium condensate 

We are now in the position to present and discuss numerical results for the in-medium quark 

condensate. Fig. 5 shows the condensate ratio (^q')(p)/(0|gg|0) as a function of the nucleon 
density p = 2A;|/37r^ in the region < p < 0.36 fm^'^ = 2.25po (i-e. kf < 345 MeV) for three 
different values of the pion mass, m.,^ = (0, 70, 135) MeV. The dashed line in Fig. 5 corresponds 
to the linear density approximation using the empirical central value cr at = 45 MeV of the 
nucleon sigma-term. One observes a very strong and nonlinear dependence of the "dropping" 



^In principle, this scale dependence is balanced by the parameter C(A) in eq.(27). But this (formal) point of 
view introduces the need to fix the scale A in an estimate of C. 




Fig. 6: Interaction contributions to the ratio between in-medium and vacuum chiral conden- 
sate. The five classes described in subsections 2.1-2.5 are consecutively added in the sequence: 
linear ^ Itt — >• iterated ^ 27r — >• A — >• ci. 



condensate on the actual value of the pion mass m^. In the chiral limit, m.,^ = 0, the quark 
condensate decreases effectively with a slope 1.8 times as large as given by the linear density 
approximation. As a consequence chiral symmetry would be restored already at about 1.5po if 
the up- and down-quark masses were strictly zero. This faster decrease is caused by two features: 
First, the ratio (jjv/m^ = — 4ci = 3.72 GeV~^ is, in the chiral limit, about 1.5 times larger than 
at the physical point, 45/1352MeV-^ = 2.47 GeV-^ Secondly, the two -pion exchange effects in 
the chiral limit drive the condensate ratio further down. Of course, the actual density at which 
chiral symmetry restoration would occur can only be roughly estimated in our calculation: 
once the chiral condensate becomes too small the very foundation of the chiral effective field 
theory approach to nuclear matter (namely the spontaneous breaking of chiral symmetry in the 
vacuum) is lost. 

At the physical value of the pion mass, = 135 MeV, the density dependence of the con- 
densate ratio (qq) (p) / {0\qq\0) is drastically different. At densities around 1.8po the in-medium 

condensate stabilizes now at about 60% of its vacuum value, and there is no further reduction in 
the entirely density region where the present chiral approach to nuclear matter can be trusted. 
For higher values of the pion mass the effects counteracting chiral restoration become even more 
pronounced. 

Let us have a closer look at individual contributions. At 2po = 0.32 fm~^ (corresponding 
to kf = 331.4 McV) one gets from the sequence of the five classes of interaction contributions 
(described in subsections 2.1 to 2.5) a total correction to (qq) (p) / {0\qq\0) beyond the linear 
density approximation of (0.14 — 0.83 + 0.27 + 1.34 — 0.54) = 0.38. Notice the cancellation 
between large contributions of opposite signs from iterated l7r-exchange and 27r-exchange with 
virtual A(1232)-excitation. At normal nuclear matter density po — 0.16 fm~^ the individual 
entries are about a factor 4 smaller: (0.04 — 0.24 + 0.07 + 0.33 — 0.13) = 0.07, suggesting an 
approximate p^-dependence of the total interaction contribution. 



Fig. 6 shows separately the effects of the five classes of interaction contributions. They are 
consecutively added up in the sequence: linear — > Itt — > iterated — > 27r — > A — > ci. The last 
two major steps (in opposite directions) should not be misinterpreted as a sign of bad conver- 
gence since the corresponding terms belong to the same order in the small momentum expansion 
of the nuclear matter energy density (— 4ci ^ 3g\/4A). Fig. 6 actually demonstrates the promi- 
nent importance of the 27r-exchange interaction beyond leading order for the in-medium quark 
condensate. 

It is important to include all effects generated by the pion-loops. For example, if one would 
drop the last constant term proportional to 4m^ ~ m^/^ in eq.(15) an amount of 0.25 would be 
missing in the total balance (at po)- As a consequence of that omission the condensate ratio 
would lie appreciably below the linear density approximation. This particular non-analytical 
4m^-term gives also an explanation for the drastically different behavior of the in-medium 
condensate in the chiral limit and for the physical pion mass. Its contribution to the condensate 
ratio, 27g'^m^p^/(1287r/^A), vanishes in the chiral limit, = 0, but reaches 100% for the 
physical pion mass = 135 MeV at p = 0.32 fm~^ = 2po- This selective consideration does, of 
course, not mean that the other one- and two-pion exchange contributions would not also depend 
strongly on the pion mass. Their dependences can be studied easily case by case with the help 
of the analytical formulas presented in section 2. In general, one can say that the condensate 
ratio {qq)ip)/ {0\qq\0) is affected significantly by interaction terms which otherwise play only a 
marginal role for the nuclear matter equation of state E{kf) (as e.g. the chiral vrvrA^A^-contact 
term proportional to Ci). This pronounced shifting of weights comes from taking the derivative 
with respect to the squared pion mass m^. On the other hand, the Fermi gas approximation (i.e. 
the linear density approximation) works reasonably well for the in-medium chiral condensate 
almost up to nuclear matter saturation density po — 0.16 fm"^, even though it is far from 
correctly describing nuclear matter as a self-bound Fermi liquid. 

Let us compare our results for the in-medium chiral condensate with other works which 
have treated to some limited extent the pion-exchange dynamics in nuclear matter. In the work 
of Lutz et al. [8] lyr-exchange plus an adjustable NN-contact interaction have been iterated to 
second order. In the region < p < 2po they find a weaker (positive) deviation from the 
linear density approximation, with no trend for stabilization of the in-medium condensate. This 
difference comes from not including the irreducible 27r-exchange, the chiral TrA^A-dynamics and 
the Ci-contact vertex. Recently, the Tiibingen group [18] has employed the chiral nucleon- 
nucleon potential at next-to-leading order to calculate nuclear matter within the relativistic 
Dirac-Brueckner-Hartree-Fock approach. They also derive the in-medium chiral condensate 
by exploiting the Feynman-Hellmann theorem cq.(l). Irrespective of using the Hartrcc-Fock 
or Brueckner-Hartree-Fock approximation, their effects from pion-exchange interactions (e.g. 
~ 15% at 2po [18]) are much smaller than in our calculation. Moreover, there is no trend of 
stabilization in the whole density region < p < 3po considered. Again, these substantial 
differences arise in ref.[18] from taking only the next-to-leading order chiral NN-potential (i.e. 
Ivr-exchange and irreducible 27r-exchange) , but neglecting the actually more important effects 
from 27r-exchangc with virtual A-isobar excitation. Wc note that our pcrturbativc calculation, 
when truncated to l7r-exchange, iterated l7r-exchange and irreducible 27r-exchange exchange, 
would lead to an in-medium condensate below the linear density approximation (see Fig. 6) . 
Qualitative differences may also come from the approximations (angle-averaged Pauli- blocking 
operator, etc.) used in the Brueckner-Hartree-Fock calculation of ref.[18] and not treating the 
iterated l7r-cxchange in full detail (as done in the present work). Apart from all these differences, 
ref.[18] agrees with our conclusion that the short-range NN-dynamics plays only a minor role 
for the in-medium chiral condensate. It is indicated to repeat the calculation of ref.[18] with the 
more complete next-to-next-to-leading order chiral NN-potential including the non-analytical 
polynomial pieces. 



_l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

^0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Fig. 7: The effective in-medium nucleon sigma-term aN,es{p) versus the nucleon density p with 
its error band ±8 Me V. 



Our results for the in-medium chiral condensate can be alternatively summarized by the 
density dependent effective nucleon sigma-term: 

a.M = - Io4 + 664) ^ ■ 

It captures the correlation effects in the nuclear medium which reduce the tendency towards 
chiral symmetry restoration. Fig. 7 shows the effective nucleon sigma-term^ as a function of the 
density for < p < 0.36 fm~"^. In this figure we display also the uncertainties associated with 
the still existing error band of ±8 MeV in the empirical determination of a^. Moreover, we have 
varied the regularization scale A between 0.6 GeV and 1.2 GeV which leads to some widening of 
the error band as the density increases. 



4 Summary and concluding remarks 

In this work we have used in-mcdium chiral perturbation theory to calculate the quark con- 
densate {qq) (p) beyond the linear density approximation. The pertinent correction term follows 
from differentiating the interaction contributions to the energy per particle of isospin- symmetric 
nuclear matter with respect to the pion mass. Analytical expressions for the contributions to 
D{kf) — dE{kf) / drn^ from Ivr-exchange (with m^-dependent vertex corrections), iterated Ivr- 
exchange and irreducible 27r-exchange with inclusion of A-isobar excitations and Pauli-blocking 
corrections have been presented in section 2. 

We find a strong, nonlinear dependence of the " dropping" in-medium condensate on the value 
of the pion (or light quark) mass. In the chiral limit, m^r = 0, chiral restoration seems to be 
reached already at about 1.5 times normal nuclear matter density. By contrast, for the physical 



qualitatively similar result, though based on a different approach, has been reported in ref.[19]. 



pion mass m^r = 135 MeV, the in-medium condensate stabilizes at about 60% of its vacuum value 
above that same density. Including systematically the effects from 27r-exchange with A(1232)- 
isobar excitation (or the equivalent chiral inrNN contact interactions 02,3,4) is crucial in order to 
obtain such a pronounced behavior. Non-analytical (contact) terms in the quark mass (such as 
~ m^^"^), which are often dropped in the presentation of the chiral NN-potential, have a very 
strong influence on the in-medium chiral condensate. Below 3po/4 = 0.12 fm~^ the correction 
beyond the linear density approximation remain relatively small. This finding can be taken 
as an a posteriori justification of the assumptions made in ref. [20] about the in-medium scalar 
mean-field. 

As a consequence of the hindered tendency towards chiral symmetry restoration (in the real 
world with = 135 MeV), pions and nucleons can be used as effective low-energy degrees of 
freedom at least up to twice nuclear matter density. The major source of uncertainty for the 
density dependence of {qq){p) is caused by the error band in the empirical determination of 
the nucleon sigma-term aN — (45 ± 8) MeV. One can hope that upcoming dispersion relation 
analyses of vrA^-scattering data and lattice QCD calculations will lead to a more accurate value 
of (T^r. Of course, there remain also questions about the size of effects from yet higher order 
interaction contributions related to 37r-exchange, two-loop 27r-exchange etc. On the other hand, 
an estimate based on recent lattice QCD results indicates that the short-distance NN-dynamics 
has only a very small effect on the density dependence of the quark condensate (qq) (p) . 
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